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1. Introduction 

The box-ball system (BBS), proposed by Takahashi and Satsuma [9], is one of the 
most important cellular automata obtained from discrete integrable systems through 
a limiting procedure called ultradiscretization [T2]. It is well known that the time 
evolution of the original BBS is determined by the ultradiscrete KdV (u-KdV) 
equation: 

= min ( 1 - [/W, J2 (L7f - uf ) , n,t € Z, (1) 

where <E {0,1} denotes the number of balls in the nth box at time t. The 

equation ((T|) is an ultradiscrete analogue of the discrete KdV (d-KdV) equation, and 
consequently the original BBS has ultradiscrete soliton solutions. It is also known that 
the nonautonomous discrete KP (nd-KP) equation with 2-reduction condition yields 
a time evolution equation of the BBS with three extensions pQ: box capacity |10j . 
carrier capacity [5] , and kind of balls [IT] . 

In this paper, we consider another type of time evolution equations of the BBS 
first presented by Nagai et al [5]- They discovered that the ultradiscrete Toda (u- 
Toda) equation on a (non-periodic) finite lattice 

( n n—l , 

^ = K^EQf-EQrM, (2a) 

j=0 j=Q 
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Et»V=E$)-Q<£?+QW, (2b) 

= E% ] = +oo, (2c) 

which we refer to as the finite u-Toda lattice simply, determines the time evolution 
of the original BBS. In this equation, the dependent variables Qn and denote 
the size of the nth soliton at time t and the size of the nth empty block at time t, 
respectively. Figure [T] shows an example of the connection between the finite u-Toda 
lattice and the BBS. Wc call this representation finite Toda representation of the BBS. 



L i Ui L 2 rjt'V'WV'V) 1 '' 

I II II II II 1 -"l «I n % «2 

t=0: .11111 1111. ..11 5 5 4 3 2 

1: 11111. ...111. .111 5 4 3 2 3 

2: 1111. ..11. ..11111 4 3 2 3 5 

3: 111. .111 11111 3 2 3 5 5 

4: 11. ..1111 11111 2 3 4 6 5 

5: 11 1111 11111 2 5 4 7 5 

6: 11 1111 11111.. 2 7 4 8 5 



Figure 1. Example of the connection between the finite u-Toda lattice and the 
original BBS. The variables Q„ and E% denote the size of the nth soliton and 
of the nth empty block at time t, respectively. 



The correspondence between the u-KdV equation and the finite u-Toda equation 
via the BBS is similar to the Euler-Lagrange correspondence of cellular automaton [3]. 
This terminology comes from hydrodynamics; the dependent variables of the Euler 
representation denote the number of particles at each point and the ones of the 
Lagrange representation denote the position of each particle. According to these 
definitions, we use the following terms in this paper: 

• Euler representation of BBS: the equation of the BBS with the variables which 
denote the number of balls in each box. 

• Lagrange representation of BBS: the equation of the BBS with the variables which 
denote the start position of each soliton and each empty block. 

• Finite Toda representation of BBS: the equation of the BBS with the variables 
which denote the size of each soliton and each empty block. 

The u-KdV equation (flj and the finite u-Toda lattice (J2J) gives the Euler representation 
and the Lagrange representation of the original BBS, respectively. Additionally, if 
we know the start position of the first soliton, we can calculate the start positions 
of all solitons and empty blocks from the values of the variables of the finite Toda 
representation. In other words, the finite Toda representation and the Lagrange 
representation can be transformed to each other. 

It has not been clarified by now why these two different ultradiscrete equations 
with different type boundary conditions describe the same original BBS. Moreover, the 
finite Toda representation of the extended BBSs are not studied sufficiently. Among 
the three extensions for the BBS, two types of the finite Toda representations have 
been clarified: Tokihiro et al [11] discussed the case of with several kinds of balls using 
the finite ultradiscrete hungry Toda lattice; Tsujimoto and the author [5] showed that 
the finite nonautonomous ultradiscrete Toda (nu-Toda) lattice determines the time 
evolution of the BBS with a carrier. The main purpose of this paper is to discuss the 
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remaining case; we derive a finite Toda representation of the BBS with box capacity. 
For this purpose, we use a map from a state of the BBS to a binary sequence, which 
we call "expansion map" . By using the expansion map, we can define the size of 
solitons in the BBS with box capacity for all time t, especially for interacting solitons. 
We also consider the finite Toda representation of the BBS with both box capacity 
and carrier capacity using new carrier rules "size limit for solitons" and "recovery of 
balls" . Furthermore, we present a particular solution for the fixed box capacity case. 

The outline of this paper is the following. In section [2j we recall the derivation of 
the Euler representation of the BBS with a carrier from the 2-reduced nd-KP equation. 
In addition, we introduce new carrier rules "size limit for solitons" and "recovery of 
balls" which are used to derive the finite Toda representation of the BBS with a 
carrier [5] and play an important role in section SI In section [3J we recall the finite 
Toda representation of the original BBS and give some remarks. In section |4l we 
discuss the finite Toda representation of the BBS with variable box capacity. First we 
discuss the case of variable box capacity A„ and no restricted carrier capacity. After 
that, we discuss the case of both box capacity and carrier capacity are variable. In 
section [SJ we give a particular solution to the finite Toda representation of the BBS 
with fixed box capacity. In section [6j we give concluding remarks. 



2. Euler representation of the BBS with a carrier 



The nd-KP equation is given by [T3] 

(an-fitJ/^V^^+^-cO/^/^'^+^-anJ/^+V^i'* = 0, k,n,t€ Z.(3) 
It is shown that the iV-soliton solution to the nd-KP equation (|3) is presented by 



J n 



fc,t 
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JC{0,1,...,JV-1} 



n w ^ n h 



i i,jeJ 



ieJ 



fe.t 
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fc-1 
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l ' 3 ' (Pi ~ Qj){9i -PjV 



3=0 J " j=0 

where ^, pi and Qi, i = 0, 1, . . . , N — 1, are some constants. Now we impose the 2- 
reduction condition with respect to the variable k, that is /^ +2,t = and Ck+2 = Ck 
for all k S Z, and set a n = 1 + S n , bt = —fJ>t, cq = 1 and c\ = 0. Then the nd-KP 
equation (|3]) reduces to the forms 

(1 + S n + pt) fig 1 /? = (1 + t*)fntJn' t+1 + Snf^fl'tl, (4a) 

(1 + 5 n + M'fZli 1 = (1 + Sjfi&W + »tft t+1 fn'ti, (4b) 
and an ./V-soliton solution to the reduced equations is given by 



fn' = 1 



( ^ 



JC{0,1,...,JV-1} 



i,j&J i&J 



k = 0,l, 



(5a) 



/ 



3=0 



Pi 



-n 1 + Mi - Pi 



1 — Pi ,0,t 



J=0 
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Wi 



u i,3 



1 -Pi~ P 3 



Let us define the dependent variables as 



It) _ Jn+l 



0,4+1 f l,t+l 



Jn+l Jn 



f O,tfO,t+l 
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Then the 2-reduced nd-KP equation (Q} yields the equations 
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(5c) 

(6) 

(7a) 
(7b) 

(7c) 



z n ' = c _z " f) / e , S n = e _A '*/ e , fi t = e~ M *A into © and taking a limit e 
obtain the ultradiscrete system 



holds. For positivity, we choose the parameters as < 6 n < 1 and < fit < 1 for all 
n, t £ Z. When the values of the dependent variables are all positive for all n, t G Z, 
we can ultradiscretize the equations J?]): putting v$ — e~ u « ) / e , v,n = e _c/ ™ A j 

fO, we 

(8a) 
(8b) 
(8c) 



tt(*+1) 
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n+1 ) 



where A ra ,M t > for all ra,tgZ. Note that we have used the fundamental formula 
for the ultradiscretization: 



lim — elog(e~ 

e->+0 



An TV-soliton solution to the ultradiscrete system ((5J is obtained as follows. Let us 
take the constants pi, i = 0, 1, . . . , N — 1, to satisfy the condition < pi < 1. Putting 
-" - ,r ' ' & = e-^A, Wij = e-^/ £ into © and 

([B]), and taking a limit e - 



, Pi = e 
0, we obtain 
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and Pi > 0, i = 0,1,..., N- 1. 

Let us introduce the time evolution rule of the BBS with the nth box capacity 
A n and the carrier capacity M t +i from time t to t+ 1. We consider the time evolution 
rule from time t to £ + 1 as the composition of size limit process and recovery process. 

(i) Size limit process: the carrier of balls moves from left (n = — oo) to right 
(n = +oo). When the carrier passes each box, the carrier gets all balls in the box; 
and if the number of balls exceeds the carrier capacity Mt+i, the excess balls are 
removed from the system. At the same time, the carrier puts the balls which the 
carrier holds into the box as many as possible. 

(ii) Recovery process: after the size limit process, all the removed balls are recovered 
to the boxes in which the balls were. 

Figure [2] illustrates these rules. 



Removed 




Carrier 



nth box 



f('+i) 



u 



(t+1) 



Figure 2. Illustration of the time evolution rule of the BBS with a carrier. 
The left figure illustrates the size limit process ]8aH and llSbl l. and the right one 
illustrates the recovery process IIScl l, 



Suppose that the dependent variables Un \ U^ +1 ^ and Z^ +1 ^ denote the following 
quantities: 

• Un^ G {0, 1, ... , A„}: the number of balls in the nth box at time t; 

• U^ +1 ^ G {0, 1, ... , A„}: the number of balls in the nth box after the size limit 
process from time t to t + 1; 

• G {0, 1, . . • , M t+ i}: the number of balls in the carrier arriving at the nth 
box in the size limit process from time t to t + 1. 

Then the equations (JSJ) give the time evolution rule: the equations (|5a|) and (|8b|) define 
the size limit process 

m (t) } + °° h> at7 (t+1) } +00 sz (t+1) \+ x ) 

I ^ n in— — oo ' VL^n Jn--co' I n J n— — ooJ 
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U ( n +1) + ((Z ( n +1) + UP) mm(zi t+1) + U^\M t+l) ) 
= ^ t+1) +* +1) + ^ ) -^ + + 1 1) -^ +1) 

gives the number of removed balls by the size limit at the nth box, the equation ()8c[) 
defines the recovery process 

,rjj{t+l), +00 r-^(*+l)| + oo \ r r/ (t+i)-i+oo 
Vl^n in— — oo> I n in— — ooJ \ w Tl in— — oo • 
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Figure 3. Example of a 3-soliton solution to the time evolution equation j5) . 
The leftmost box is the Oth box in the both figures. The carrier capacity Mt = 6 
for all t > 1 . Each number denotes the number of balls in a box and ' . ' denotes 
an empty box. In the right figure, boxes containing recovered balls are shown in 
boldface (compare to the left figure). 



Figure [3] shows an example of a 3-soliton solution to the time evolution equation 
([5]). The parameters are chosen as 

{3 if n is even, f +oo if t < 0, 

M t = I 
5 if n is odd, 1 6 if t > 0. 

Remark 1. Eliminating the variable U^ +1 ^ from the equations (|8a[) and (|8b[) . we have 
the relation 

- = min(A„ - [/W, Zi t+1) ) - min(zi t+1) + U^,M t+1 ). (9) 

From the equation (|5c| . the relation 

n-l 

z { r ] = E (c/f -c/j t+1) ) (io) 

J=— oo 

also holds. Substituting (fTQ)) into (|9]), we obtain the equation 

n-l 

r(*+i) - ™^ I a _ rr(*) V" - n^h 
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v , J= ~~ ^ , (in 
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+ max 0, U f ] ~ E ^ t+1) -^+i 



J — — OO J — — OO 
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where we have used the formula 

- min(-yl, —B) = max(A, B). 

The equation (TTT]) has the same form as of the time evolution equation of the "BBS 
with a carrier" presented by Takahashi and Matsukidaira [5]. 
If we choose M t +\ = +00, then (fTTj) yields 

Uf l) = min L - J2 (uf - U^)) , (12) 

\ 3'=-OQ J 

which is the nonautonomous u-KdV equation. In addition, (|8b[) yields 



_ ^(t+i) (t) -(t+i) 



(t+1) _ tt(*+1) 



and comparing this relation with (I5cl) , we have the relation L^A = U 
3. Finite Toda representation of the original BBS 

We recall the relation between the finite u-Toda lattice and the original BBS. The 
bilinear form of the discrete Toda lattice is given by 

^-M +1) =^- 1 1) e + t 1) +« ) , M€Z. (13) 
Let us introduce the dependent variables 

T W T (t+i) It) (t+i) (t) (t) 

(t) _ r " T »+l (t) _ T n+l T rt-l ,(t) _ T " T 7t+1 

g " J*) C « JJt+l) ' Jt-l)Jt+l)' 

'n+l'n 'n+1 '« 

Then (JT3J) yields the equation 

9i t+1) =ei*l 1 + ^ +1 ), (14a) 
and the identities 

(*) (*) 

c n c n (t+1) ' ™ "n— 1 (t+1) \ LW I 

1n-l Qn-1 

hold. Putting q® = e^'A e W = e^'A ^ = e- ^, and taking a limit 
e — > +0, we obtain the u-Toda lattice 

Q(*+ 1 )=min(^j 1; ^+ 1 )), (15a) 

E^=eU-Q^+QU, (15b) 

^ t+1) = ^n-l 5 - Qn-1 1 + Q [ n- (I5c) 

Furthermore, we impose the terminating condition for discussing the finite Toda 
representation: 

E®= ^=+00, d£ +1 >=qW, (15d) 

where iV is a positive integer, which denotes the number of solitons in the original 
BBS. 

Let the variables Qn\ E%' and Dn +1 \ respectively, denote the following 
quantities of the original BBS: 



Qn : the size of the nth soliton at time t (n = 0, 1, . . . , N — 1); 
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• E$: the size of the nth empty block, namely, the distance between the (n — l)th 
soliton and the nth one at time t (n = 1, 2, . . . , N — 1); 

• Dn +1 ^'- the number of balls in the carrier after getting Qn^ balls (n = 0, 1, . . . , N— 
!)■ 

Then the next theorem gives a fundamental result on the connection between the finite 
u-Toda lattice and the BBS. 

Theorem 1 (Nagai et al 5 ). The finite u-Toda lattice H5\) determines the time 
evolution of the original BBS. 

Remark 2. Conventionally, the d-Toda lattice (qd-type) has been written in the form 

„{t+i) , Jt+i) _ (t) , At) (t+i) (t+i) _ (t) (t) 

'in ~ ~ In Tc n+li Hn—1 n ~ In c n J 

or 

g (W)= g «-e(f+ 1 )+ei t | 1) (16a) 
ei t+1) =e«^, (16b) 

9n-l 

which we cannot ultradiscrctize directly due to "negative problem" . The equations 
(|16l) are called Rutishauser's qd algorithm in numerical algorithms [B]. On the 
other hand, the equations (fH]l are called dqd algorithm. The dqd algorithm is the 
subtraction-free form of the qd algorithm and computes matrix eigenvalues or singular 
values more accurately than the qd algorithm. 

Suppose the finite lattice condition eg = e$ = 0. Nagai et al [5] rewrote (|16ap 
using (|16b|) as 

o (t+1) -e (t) + o (t) -e (t+1) 

Hn ~ e n+l • In e n 
(*) 

- Jfi) g " (g {t+l) c (t) ) 

<?n-l 

- e (*) | g " } f («) _ e (t+1) ) 



e 



(i) . llj=ogj 



™+l ^ TTfl-1 (t+1) ' 

n J= o q 

Then they could ultradiscretize the finite Toda lattice: 

(n n — 1 \ 

^Esf-E^h ( 17a ) 

^+ 1 )= J EW-Q^+QW, (17b) 

= = +oo. (17c) 

On the other hand, by introducing an auxiliary variable 
At+i) ._ (t) _ (t+i) _ (t+i) _ (t) 

we can ultradiscretize the finite d-Toda lattice directly without the negative problem 
and obtain the finite u-Toda lattice of the dqd form ([151) . From the viewpoint of the 
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BBS, the variable Dn denotes the number of balls in the carrier. Therefore, the 
finite u-Toda lattice of the dqd form (fT5j) is important to consider the finite Toda 
representation of the BBS with a carrier. 

Remark 3. Here we remark on the Lagrange representation of the BBS, which is also a 
terminology from hydrodynamics as the Euler representation; the dependent variables 
of the Lagrange representation denote the position of solitons. Let the variables X„ 
and Yn^ denote the start position of the nth soliton and the one of the nth empty 
block at time t, respectively (see Figure H]). Then the Lagrange representation of the 
BBS is given by [I] 
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n+1' 
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Fj^) = yW + min ( X® £(lf - Xfl x ) - J2(y} t+1) ~ X<F») ] , (18b) 
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Relations between these variables and the variables of the finite u-Toda lattice (fTTf 
are given by 

x« = Y« + i#), yw = x& + (is) 

We can readily show that (fT8|) and (fT9| yield the finite u-Toda lattice (ITTl) . Conversely, 
we can calculate the values of {X}i }n=i an d {y£ }n=i Ir0m the values of Xq, 
{Qn }n=o an£ l {-^" }n=i usm g the relations P^|. In other words, the finite u-Toda 



lattice (fl7|) and the equation Xq 



(*+i) 



(|19p . uniquely determine the time evolution of the BBS. 



• Qq , which is obtained from (|18a|) and 
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Figure 4. Lagrange representation and finite Toda representation of the BBS. 



4. Extension of the finite Toda representation to the case of variable box 
capacity A„ and variable carrier capacity M t 

In previous studies, the finite Toda representation is considered only for the BBS with 
box capacity 1. In this section, we extend the finite Toda representation to the case 
in which each box has own capacity A„. First we consider the case of carrier capacity 
M t = +oo. The Euler representation of this case is given by (|12p . 

Wc first define the size of solitons and the one of empty blocks for the BBS with 
variable box capacity A„ at any time t. For this purpose, we refer to the work of 
Takahashi and Satsuma [TU]. They analyzed the BBS with the fixed box capacity A 
using a map from a state of box capacity A to a binary sequence. We generalize this 
map for the case of variable box capacity A n . 
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Suppose that a state of the Euler representation (fT2|) {Un^}n=-oc such that 
UrP G {0, 1, ... , A„} is given. We assume that, for simplicity, Un — for n < 0. Let 



us define a map {ufi}*™^ i-> {V^}^^, where V^' £ {0, 1}, as follows: 

(1) K (t) = for n< 0. 

(2) Let s = and s„ = X^Co f° r n = l,2, From n = to +oo, if V^Li = 1, 

then 
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Note that the relation U, 
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s „+A„-i nQlds _ 



We refer the jth number vj in the binary sequence as the jth segment. By 
using this map, the nth box is expanded to the block composed from the s„th to 
(s„ + A„ — l)th segments. We call this map expansion map from a state of the BBS 
with variable box capacity A„ to a binary sequence. 



.1111111111. 



A: 5353535353 
U: .15213. .2. 

I 

...111111 I 



.11 



8 5 4 11 2 



Figure 5. Example of the expansion map from a state of the BBS with variable 
box capacity A n to a binary sequence. In the binary sequence, a block between 
two separators ' I ' corresponds to an original box. 
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Figure 6. Example of the expansion map for the BBS with variable box capacity 
A n . The carrier capacity parameter is chosen as Mt = +oo for all t £ Z. 



Figures [5] and [S] show examples of the expansion map. As shown in Figure [HI the 
expansion map enables us to define the size of the nth soliton Q„ and the one of the 
nth empty block En for the BBS with box capacity A„ at any time t in the same way 
as for the BBS with box capacity 1. Let D% +1 ' denote the number of balls which the 
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carrier holds after getting Q$ balls, and A$ denote the capacity of the box which 
contains the beginning (leftmost) segment of the nth empty block. Then we arrive at 
the following theorem. 

Theorem 2. Let the variables Qn\ Er? and D n t+1 ^ denote the quantities of the BBS 
as explained in the previous section. Then the time evolution of the BBS with variable 
box capacity A« is given by 

(20a) 
(20b) 



= min (E^ +1 max(0, ^ - D^), 

c(0, A« - d£?) + max(0, A® - #£ t+1) ), 



— max 



£)(«+!) 



E, 



(*) 



E 



D 

(t) 



(t+i) 

n-1 



-oo, Dt 1] 



)(*) 

5n ' 



(20c) 
(20d) 



We note that, from (|20a|) . Qn < Dn holds for all n and t. Since the size 
of the nth soliton Qn ^ should be positive for all n and t, from (|20c[) , the inequality 



Dn +1> > 1 holds for all n and t. Thus, all max(0, A^ 1 ^ — D n l+1> ) are equal to zero 

when An — 1 for all n and i, the case of the original BBS. In this case, the equations 
(|20p reduce to the finite Toda representation of the original BBS (TT5)) . Hence we can 
say that the ultradiscrete system (|2"0"|) is a generalization of the finite u-Toda lattice 

IIS]). 



(0 



Proof. Let us consider the general An^ > 1 case. As we mentioned above, (|2T)|) has 
additional terms max(0, A^ +1 — D n t+1 ^ ) which do not appear in the case of box capacity 
1 (|T5|) . Hence, we shall investigate the role of the terms max(0, A^ +1 — D„ +1 ^). 



a-" 



- H - | 



max(0,A™, --d;; +1) ) ■ 



■ £/<,' +11 = min(Aj£, - Bj>, £><' + " - £/<,") 



Figure 7. Illustration of the quantity max(0, A^ij — -D£* +1 '). Wc now need to 



determine the quantities on the area filled with gray: K 



(*+i) nC+l) 



In this figure, - 
just before getting Q„ balls. 



1 is assumed; the carrier has one ball 



Let us consider the time evolution of the BBS with box capacity A„ from time 
ttat+1. Assume that Q { - +1 \ j = 0, 1, . . . , n - 1, and Ef +1 \ j = 1,2,..., n- 1, 
are given (see Figure [7]) . Let m be the index of the box which contains the leftmost 
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segment of the (n + l)th empty block at time t. Then the capacity of the mth box 
A m is equal to A^*^ by definition. Moreover, the relation 

n Tl—l m— 1 

= £ Qf - £ Qf +1) = £ (fjf - c/j t+1) ) + u<$, 

where ujp denotes the number of balls in the kth box at time t, also holds by definition. 
Hence, we can calculate the quantity Um by the nu-KdV equation ([1 

= min L m U$ t £ (C/f - uf^) 

\ j=-oo 

= min(Ai t | 1 -^ ) , J Di t+1) -^ ) ). 
Then we obtain the relation 

A - uW - r/( t+1 ) - A (t) - r/ (4) - miniA (t) - f/ (t) n( t+1 ) - r/Wi 

= -min(0,Dl*+ 1 )-Ai*{ 1 ) 

= max(0,Ai t | 1 - J D(*+ 1 )), 

where we have used the identity — min(— A, —B) — max(A, B). This relation implies 
that the term max(0, J^}+i ~ Dn +1 ^) denotes the size of interspace inserted between 
the nth soliton at time t and the nth one at time t + 1 . 

Once we notice the role of the terms max(0, A^^ — Dn +1 ^), we can now clarify the 

meaning of the equations (|20aD and (|20bp . Since the term — max(0, A^ +1 —Dn +1 ^) 



in (|20ap denotes the difference between the size of the inserted space and the one of 
the (n + l)th empty block, Qn should be determined by (|20a[) . Similarly, En +1 ^ 
should be determined by (|20bl) . It is also true for n = and 1, then the proof is 
completed by induction. □ 

Next, we construct the finite Toda representation of the BBS with both box 
capacity and carrier capacity from two time evolution maps: the size limit map and 
the recovery map. This is a similar as for the construction of the Euler representation 
explained in section [2l Figure [8] shows an example. 

In the next theorem, we use the following notations: 

. Q { n\ E$: the size of the nth soliton and the one of the nth empty block at time 
t, respectively; 

• Qn +1 \ : the size of the nth soliton and the one of the nth empty block 
after the size limit process from time t to t + 1; 

• C„ +1 \ D^ +1 ^ : some quantities which will be explained in the proof of the next 
theorem in detail; 

• Kn\ a!*': the capacity of the box which contains the leftmost segment of the 
nth soliton and the one of the nth empty block at time i, respectively. 
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Figure 8. Example of the expansion map for the BBS with box capacity A n and 
carrier capacity Mt = 6 for t > 0. This is obtained from the example in Figure [3] 
via the expansion map. 



Theorem 3. Let the variables Qn\ En \ Q„ + ^ > E^ +1 \ C^ +1 and D^ +1 denote 
the quantities of the BBS as explained in the above. Then the time evolution of the 
BBS with box capacity Kn^ and An\ and carrier capacity M t +i is given by 

Q ( : +1) = min (£gx - max(0, -D^D^) , (21a) 

^ +i) =^ ) -e + 1 i) +^ (21b) 

- max(0, A# - Z& 13 ) + max(0, Ag. x - 

Ci t+1) = min(Si'+ 1) - + K%\M t+1 ), (21c) 

Z>i t+1) = min(ci t+1) + Q<*> - K£\M t+ i), (21d) 

Ql t+1) = Q« + - - jtf> + 41, (21e) 

+ -Q« (21f) 

^ = = ^ = ^ = = (21g) 

where the carrier capacity M t+ \ must be chosen to satisfy the condition Kn < M t +i 
for all n and t. 

When the quantities {Q^j^o and { E n } }n=i 

are given, first we can calculate 

Dq +1) using (|21d)) and ( gig] ). Next, we can calculate by ([2"Taj) . cf +1) by (j2"Tc|) , 

!D X by If21d| . In a repetitive manner, we can obtain the quantities {Q n }n=oi 

{^n +1) }n=o and {^ ( * +1) «}^=o 1 - Finally, we can calculate the quantities {E { * +1) }^, 

{Qn +1) }n=o and {En +1) }n=i h Y WM > <EHl and (EH, respectively. Hence the time 
evolution is determined by (f2"Tj) . 

If M t+ i = +oo, then dHc|> a nd (EcQ) reduce to ci t+1) = - + K$ 

and D^ +1 = C^* +1 ^ + Qn) — Kn\ respectively. Thus we have the equation 

D ( , t l +1) = D^l^ - Qn-i + Qn } and > substituting them into (|2Tc|) and JUJ, we 

obtain = Qn +1 ^ and En +1 ^ = E* +1 . Hence, in this case, the ultradiscrete 

system (f2"Tj) reduces to the system (|2"0)) . We can therefore say that the system (|2"Tj) is 
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a generalization of the system (QUI) . 

Proof. Let us show that the equations (|21al) - (|21d|) describe the size limit process and 
(|21e[) - (|21fj) describe the recovery process. 

First, we consider the size limit process. The equations (I21a|) and (I21bp have 
the same forms as of (|20a[) and (|20b[) . Thus we shall investigate the variables 

C n * +1 '' and D^ +1 ^ which are defined by (I21c[) and (|21dl) . Suppose that the carrier 

capacity is chosen as K$ < M t +\ < +00, Q* +1 \ j = 0, 1, . . . ,n — 1, and E* +1 , 

j = 1, 2, . . . , n — 1, are given, and the quantity D^-i denotes the number of balls 

which the carrier holds after getting Q^-i balls from boxes and restricting the number 

of the holding balls to Mt+x balls. Since the inequality Qn—i^ < -^n-i holds from 
(|21a[) . it is sufficient to consider the following two cases: the case of which the carrier 

drops off all balls temporarily (D^x — Qn-i — 0) an d the case of which the carrier 

has balls just before getting Q$ balls (D^Jx — Q^-x > 0). 

(l) If D n-l - Qn-l = °> then C n = ™(O n -l ~ Qn-l + K n , M t+ l) = K n 

holds from the assumption. We should note that, in this case, the number of 
balls which the carrier holds is zero temporarily before getting Q„ balls. Thus, 
from (|21dl) . we have D^ +1 ^ — mm(Qn\ Mt+x), which indicates that the quantity 

Dn +1 ^ is again the number of balls which the carrier holds after getting Q$ balls 
and restricting the number of the holding balls to M t +\ balls. 

(ii) The case of — Qn-l* > 0- Let m ^ e ^ ne m dex of the box which contains the 

leftmost segment of the nth soliton at time t. Under the assumption, in the terms 

of the variables of the Euler representation ([8]), D^—x — Qn-l* > implies that 

^m +1 ^ — — Um should hold (see Figure IH]). Now A m = by definition. 



t : 


• 







00 


i 














W - [T« 










+1) 






q" + " C" 


t + 1: 


O 





00 




mm 


••• 



Figure 9. Illustration of the size limit process when the carrier parameter 
M t+ i = 6 and D^-i^ ~Qn-i > 0. We can see that C^ +1 ^ denotes the number 
of balls which the carrier holds after getting Um balls (white balls) and restricting 
the number of the holding balls to Mt+i, and denotes the number of balls 

which the carrier holds after getting Q^ 1 — balls (gray balls) and restricting 
the number of the holding balls to Mt+l balls. 

Hence we can write and flUd} as C { * +1) = minp^ - (Q^-i -T7„ +1) ) + 
u£\M t+1 ) and D% +1) = mm(c£ +1) + (Q$ - U$) -U^ +1 \ M t+1 ), respectively. 
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Table 1. Change of the number of balls which the carrier holds. 

State of the carrier The number of balls which the carrier holds 



Getting balls 

Putting Q„_ 1 - U rn balls D n _ x - {Q n _ 1 - U m ) 
Getting C/^ balls - Q n _ x + iQ 



Size limit to M t+ i balls 
Putting I7^ +1) balls 

Getting Q<? - 0<? balls + Q<*> - 

Size limit to M t +i balls 



^n— 1 




u n-l 


- (e + i x 


D [t+1) 




77(«+l) 
°n 




77(*+!) 
°n 


-u [t+1) 

m 


77(*+l) 
°n 


+ q1 4) - 


^ n 





Therefore, the quantity -D n * +1 is again the number of balls which the carrier holds 
after getting Qn^ balls and restricting the number of the holding balls to M t +\ 
balls. We can summarize the change of the number of balls which the carrier 
holds as Table [U 

Thus, together with the proof of Theorem [3J it is proved that (I21a,[l — (|21d[) describe 
the size limit process by induction. 

Furthermore, since the number of balls removed by the size limit process are given 

by {D n _ 1 -Q n _ x +K^')-C n and (C„ + Qii - ') - D n , we obtain 
the equations of the recovery process 

— Vn ' °n °n+l n ri < ^n+V 

which lead to the equations (|21c[) and (|21f[) , and the proof is completed. □ 



Remark 4. Furthermore, the variables and Xq* , which denote the index of the 
leftmost segment of the Oth soliton at time t satisfy the equations 

4 t+1) =xf + QW + max(0 ! A^-Sr ) ) ) 
5. Particular solution for the fixed box capacity case 



In this section, we discuss a particular solution to the ultradiscrete system (1211) with 
a special condition: all boxes have constant capacity A. 



pO,t+l_l,t-l - X <),t-l-l,t+l , _o,t=i,t 
n 'n ~ ul n+l ' n— 1 ' 'n 'n ' 
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Let us consider the bilinear equations 

(23a) 

r==0,t J,t /r ,, \_0,t=l,t I ,, =0,t+l_l,t-l fOQM 

0,* =l,t _ =0,i 1,4 , =0,4+1 M-l / 2 o N 

'n+l'n — 'n+l'n ^ ' n 'n+1 > !^ J W 

(5 - Mt^+V^ 1 + t^t 1 / = r^I 1 r^+ 1 , (23d) 

where (5 is a constant and /it is a parameter depending on t. We introduce the 
dependent variables 

-0,4+1 1,4 j- -0,4+1-1,4 

(4) _ T n+1 T n -(4) _ T n+1 T n 

_o,t+i M ' «n j tT 0,t+l T M ' 

' n ' n+1 ' n ' n+1 

0,4 — l,t+l —0,4 — l,t+l 

(t ) ^ViVl -(«) ) -r n+l r„_ 1 

e n u 0,t— 1,4+1 ' e » U ^ U ^ =0,4=1,4+1 ' 
^~n n n 

-(4) _ r r n r n "jW _ r n+l T n 

C ™ -°-0,i+l l,t-l' fl " - g _ u -0,4+1 1,4-1 • 



'n+1 



Then (I23dp yields the relation 

1 + S-'diP = (5- fx t ] 



0,4-1-1,4+1 
-1 'n+1 ' n 

=o,t+i ■ 

1 n 'n+1 



Further, (|23al) - (|23cp yield the equations 



t +1) =^Ul + S-^ +1) ) + d { : +1 \ (24a) 
-ft+i) 

= (6 - ^ t+1 )^- + fM+1 , (24b) 

?n-l 

. \— 1-rt+ii , Mt+1 



4T J = (* - m^O" 1 ^^ + i T-T , (24c) 

1 - d Vt+i 



and the identities 



,(*) 



s d n -i 



g(t+i) _ e ffl g" 1_I ~~»- + i -> 1,1) 

^ +1) = ^|^, (24e) 

C n+1 

_(t+l)-j(*+l) 

p (t+l) _ F (t+D g »-l ( 9 Af\ 
e " ~ 6 « (t) -j(t+l) ' 

1 

hold. In addition, we impose the finite lattice condition 

e«=e«=^=4>=0. (24g) 
In the bilinear equations (1231) . this condition implies 

k,t k, t — k,t — k.t pi 

T -l — T N+1 ~ T -l — T W+1 — U ' 

We assume that the constant <5 and the parameter fi t satisfy the condition 
< (it < 5 for all t £ 1. Then, putting q n t} = e~^/ e , e„* } = e~ E ^/ e , 
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and taking a limit e — > +0, we obtain the ultradiscrete system (|2"Tj) with the condition 

K { n ] = A { n = A < M t+ i for all n,teZ. 

The following theorem is proved by using a determinant identity called the Pliicker 
relation. 

Theorem 4. A particular solution to the bilinear equations \2S\) with the semi-infinite 
lattice condition r*'* 



—k.t 

T' 



T k \ = for all k,t G Z 


is given by the 


Hankel determinants 




if n < 0, 




('. 


if n = 0, 


(25a) 


1 \tk+i+j\o<i,j<n-l 


if n > 0, 






if n < 0, 






if n = 0, 


(25b) 


I \£k+i+j\o<i,j<n-l 


if n > 0, 





where and ' are arbitrary functions satisfying the dispersion relation 

t +1) = + = (Mt - <5Ki+i + £?\ n = 0, 1, . . . . (26) 

Hereafter, we choose the arbitrary functions as 



where j»i and u>i, i — 0, 1, . . . , N — 1, are some constants. Then the dispersion relation 



(|261) is satisfied and the finite lattice condition t_\ — r Ar ' +1 = = t^ +1 = holds 
for all k, t G Z. Substituting (j2"7| to (|2"5j) and expanding the Hankel determinants 
using the Cauchy-Binet formula, we obtain 

\ - I I TT Pri - Prj 1 IT Wr< n*-=o(Pri + Mi) 



T' 



0<r o <ri <...<r„_i<iV— 1 \ \0<i<j<n- 



.! (Pr, + ^(Pr, + $) / ^ (pr, + 



2 



E 

0<ro<n <...<r n - 1 <N — l 



(Pr«+*)(Pr,+*) ii (p* +<*)*+* 



,0<i<.j'<ra-l ^ * '/ j=0 



for n = 1,2, ...,N. These expressions can be ultradiscretized directly: putting 

p„ = e"^"/ 6 , w n = e~ Wn / e , t^'* = e~ T ™'*/ £ , r^'* = e~ T »' / £ , and taking a limit 
e — > +0, we obtain the next theorem. 

Theorem 5. ^4 particular solution to the ultradiscrete system \21)) with the condition 
Kn^ — An = A < Mt+i for all n,t € Z is given by 





— 1 71+ 1 


n 


_1 j_ T 1 - 1 — T 1 

' n n 


t 

f 1* 




^ J n+1 


n 1 n 


-T 1 '* 

- 1 71+1 ! 




— 1 n+l 




mht+l _rn 1 > t 

n—1 n 


fl +2A, 




^ J ri+1 


=0,t Spl.t+l 


^ t+1 +2A, 


r<(t) 


- r '* _ 

n 


n 


' n n 


: + A, 




rriO.t 

— 1 n+l 


T°' t+1 + T M - 

n ' n 


1 n+l J 
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(n-l 
V (w n + (2(n - 1 - ») - l)P u 

-(2(n - 1) +t + fc) min(P n ,A) + ^min(P r4 ,M,-m, n = 1, 2, . . . , N, 

j=0 / 

(n-l 
V (V rj + 2(n - 1 - <)P„ 

t-1 \ 

-(2(n - 1) +t + k) min(P r! , A) +^min(P r< ,Mj)J , n = 1,2,..., AT, 

j=o / 

rpk,t rpk,t rpkft rpklt , rpk.t TrT^t ri 

J- -I — J-N+l ~ 1 -1 — 1 N+l — < 00 ) J — J ~ U ' 

where Pi and Wi, i = 0, 1, . . . , iV — 1, are some constants satisfying Po < P± < . . . < 
Pn-i- 

Remark 5. There exists a Backlund transformation from the discrete system to 
the nonautonomous discrete Toda (nd-Toda) lattice: 

n (t) d g" Jt) e « 



5(l + <5-ioi t) )(l + <5- 1 el* ) )' " ^l+J-^ilja + ^e^) 



q(t) = " ^ e (t) = 

In fact, these variables have r- function expressions 

(t) = g-1 T » T n+l -(t) (S-nA- 1 n "+ 1 

q " o,t -o,t+i ' q » ^ w =o,t =0,t+l ' 

n+l n 'n+l n 

0,t — 0,t+l — 0,t — 0,t+l 

e (t) = j r n+l r n-l -(t) = (S - lu) Tn+lTn - 1 

Si 0,t_ ,t+l ' e « [U W _0,f_0,f+l • 

T n ' n ' n n 

Since the bilinear equations 



0,t-l-0,t+l _ r/r _ \ 0,t-l_0,i+l , o,t=o,t 
r n T n —0^0 H'tPn+l 7 n-l ' T n r n J 

=0,t 0,t _ r , _ „, W 0,feO,i , „. 0,t-l-0,t+l 



c— u,t u.t /r \ u,t— u.i | u,t — i— u 

dT n r n+l = ( d - MtJV T n+1 + M* T n+l T r. 

hold (these are proved by using the Pliicker relation), we have the equations 

(*) 

q n = e n+l+ d n > d n = d n-l _( t+1 ) + <?t+l, 

where 

—o t O.t r— 1 

j(*) /r s-1 T n r n+l IH 

' n+l ' n ^ L 

Additionally, we have the identities 

Jt) d (t+1) a (t) d (t+1) a ( * +1) 

" (t+1)' Mn -(t+l)_ (t+ i) ' e « e « -(t+1) ( t ) ' 

q n-i d„ q^_/ d n _ 1 q„ ; 
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Eliminating d 



from these equations, we obtain the modified nd-Toda lattice [2] 



n 




and the finite lattice condition is given by 



6. Concluding remarks 

In this paper, we have derived the finite Toda representation of the BBS with box 
capacity by introducing the expansion map from a state of the BBS to a binary 
sequence. Furthermore, we have given a particular solution for the fixed box capacity 
case. Hence we can say that the ultradiscrete system (12"TT) is integrable if the 
parameters Kn^ and A„ & r e chosen as constants. Since there is a connection 
between the ultradiscrete system (j2"T]) and the BBS with variable box capacity 
which is integrable, we expect that the ultradiscrete system (|2"Tj) of the variable box 
capacity case is also integrable and a discrete system derived through the inverse- 
ultradiscretization has determinant solutions. This problem is left for future research. 

In the proof of Theorem[31 the variables c[ l * +1 ' ) and D^ +1 ^ have played important 
roles; these variables denote the number of balls which the carrier has. Moreover, these 
variables correspond to the variables which are introduced to remove subtractions 
in the discrete equations (see Remark [2]). This result gives us a guideline for 
ultradiscretization of Toda-type integrable systems and making connections between 
these systems and BBSs. 

In 2000, Spiridonov and Zhedanov [7] proposed a Toda-type nonautonomous 
integrable system called R\\ chain, which is derived as the compatibility conditions of 
spectral transformations for some biorthogonal rational functions. By using techniques 
developed in this paper, we will be able to ultradiscrctize the R\\ chain and consider 
a corresponding BBS. 
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